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SUMMARY 

Advances in stem cell manufacturing methods have made it possible to produce stem cell-derived cardiac myocytes at industrial scales for 
in vitro muscle physiology research purposes. Although FDA-mandated quality assurance metrics address safety issues in the manufac- 
ture of stem cell-based products, no standardized guidelines currently exist for the evaluation of stem cell-derived myocyte functionality. 
As a result, it is unclear whether the various stem cell-derived myocyte cell lines on the market perform similarly, or whether any of them 
accurately recapitulate the characteristics of native cardiac myocytes. We propose a multiparametric quality assessment rubric in which 
genetic, structural, electrophysiological, and contractile measurements are coupled with comparison against values for these measure- 
ments that are representative of the ventricular myocyte phenotype. We demonstrated this procedure using commercially available, 
mass-produced murine embryonic stem cell- and induced pluripotent stem cell-derived myocytes compared with a neonatal mouse ven- 
tricular myocyte target phenotype in coupled in vitro assays. 



INTRODUCTION 

In response to widespread efforts to commercialize differen- 
tiated stem cells (Brower, 1999), the U.S. Food and Drug 
Administration (FDA) established a set of regulations and 
guidelines for manufacturing and quality-control evalua- 
tions of human cellular and tissue-based products derived 
from stem cells (FDA, 2011). The recommendations out- 
lined for evaluating differentiated stem cell phenotypes 
were developed specifically to address patient safety 
concerns such as tumorigenicity and immunologic incom- 
patibility, due to the initial focus of the industry on regen- 
erative-medicine applications (Fink, 2009). Concerns about 
patient safety may have slowed the commercialization of 
regenerative therapies (Fox, 2011), but the use of industrial 
stem cell-based products for in vitro research, particularly 
pharmaceutical screening applications (Placzek et al., 
2009; Rubin, 2008; Thomson, 2007; Wobus and Loser, 
201 1), is a promising goal that can potentially be reached 
in the near term. Due to the mandate to test all drug com- 
pounds for potential adverse effects on the heart, in vitro 
cardiac toxicity screening is a particularly important appli- 
cation that has prompted the development of commercial 
stem cell-derived cardiac myocytes by a number of com- 
panies (Webb, 2009). In this context, the focus of quality 
assurance shifts from patient safety concerns to the devel- 
opment and adoption of measures that ensure these cells 
reliably mimic cardiac myocytes found in vivo. 

In order to develop quality assurance standards for assess- 
ing stem cell-derived myocyte differentiation, it is neces- 



sary to first establish the set of characteristics that reliably 
define cardiac myocyte identity. We reasoned that the 
most effective way to delineate these standards was to 
comprehensively evaluate the aspects of form and function 
that give rise to the contractile properties of cardiac myo- 
cytes in the healthy, postnatal heart (Sheehy et al., 2012). 
No standardized approach currently exists for evaluating 
cardiac differentiation. Basic characterization involves the 
use of one or more assays with stringencies ranging from 
the observation of spontaneous beating activity to electro- 
physiological recordings, and one of the most commonly 
used approaches is gene-expression profiling (Mummery 
et al., 2012). In addition to measuring the expression of car- 
diac biomarker genes (Bruneau, 2002; Ng et al., 2010), we 
also examined the organizational characteristics of the 
contractile myofibrils (Feinberg et al., 2012), the electrical 
activity that regulates myofibril contraction (Kleber and 
Rudy 2004), and contractile force output directly (Alford 
et al., 2010). Since human ventricular myocytes are not 
readily available, we utilized commercially available mu- 
rine embryonic stem cell (mESC)- and induced pluripotent 
stem cell (miPSC)-derived myocytes and compared them 
against ventricular myocytes freshly isolated from neonatal 
mice (neonate). Although humans and mice exhibit differ- 
ences in cardiac physiology, our goal was to determine the 
utility of comparing industrially manufactured stem cell- 
derived myocytes and isolated cardiac myocytes that 
possess the desired phenotype, using a multifactorial com- 
parison of high-level myocardial tissue architectural and 
functional characteristics. 
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RESULTS 

Gene-Expression Profiling of mESC, miPSC, and 
Neonate Engineered Tissues 

We previously reported the influence of tissue architecture 
on the contractile performance of engineered myocardium 
in vitro, so we began our characterization of commercially 
available mESC and miPSC myocytes by evaluating their 
response to geometric cues encoded in the extracellular 
matrix (ECM) (Feinberg et al., 2012) and measuring the 
expression of genes that are commonly used to delineate 
the cardiac myocyte lineage (Maltsev et al., 1994; Sartiani 
et al., 2007). Culturing the mESC (Figure lAi) and miPSC 
(Figure lAii) myocytes on a substrate coated uniformly 
with fibronectin (FN) gave rise to monolayers with an 
isotropic cellular arrangement similar to the arrangement 
observed when neonate ventricular myocytes (Figure 1 Aiii) 
were cultured in a similar manner. Moreover, mESC (Fig- 
ure SlAi available online), miPSC (Figure SlAii), and 
neonate (Figure SI Aiii) myocytes all assumed a pleomor- 
phic morphology when cultured sparsely on isotropic FN 
(Figure SIB), even though the neonate cardiac myocytes dis- 
played a smaller surface area than the mESC and miPSC 
myocytes (Figure SIC). A comparison of the expression pro- 
files for isotropic mESC- (Figure IBi) and miPSC- (Fig- 
ure IBii) derived tissues versus the neonate tissues revealed 
a number of significant differences associated with ion 
channel subunits and components of the sarcomere. In 
particular, the mESC tissues exhibited significantly higher 
expression of the L-type Ca 2+ channel subunit Cacnald 
(4.9-fold, p < 0.05), as well as the T-type subunits Cacnalg 
(9.0-fold, p < 0.05) and Cacnalh (42.2-fold, p < 0.05) versus 
neonate tissues. Isotropic mESC tissues also showed signifi- 
cantly lower expression of Irx4 (—9.1-fold, p < 0.001), Myl2 
(-3.2-fold, p < 0.05), and Myl3 (-3.8-fold, p < 0.01), which 
are commonly associated with the ventricular myocyte 
phenotype (Ngetal., 2010), and significantly higher expres- 
sion of the atrial marker genes Myl4 (40.2-fold, p < 0.001) 
and Myl7 (24.5-fold, p < 0.01) than the neonate isotropic 
tissues. In contrast, the miPSC isotropic tissues showed 
significant differences in expression for Cacnald (5.7-fold, 
p < 0.05), Cacnalh (27.9-fold, p < 0.001), Myl4 (14.1-fold, 
p < 0.05), and Myl7 (11.1, p < 0.05) versus the neonate 
isotropic tissues. These observations suggest that compared 
with the mESC-derived myocytes, the miPSC-derived myo- 
cytes exhibited an expression profile that more closely 
resembled the profile of the neonate ventricular myocytes. 

Based on previous studies, we recognized that the gene- 
expression profile of cardiac myocytes changes as a func- 
tion of the tissue architecture within which they are 
embedded (McCain et al., 2013). We engineered laminar, 
anisotropic myocardium from mESC (Figure ICi), miPSC 
(Figure ICii), and neonate cardiac myocytes by culturing 



them on microcontact-printed FN, where the cells sponta- 
neously formed cell-cell junctions and aligned with the 
geometric cues within the matrix to form a contiguous tis- 
sue of high-aspect-ratio cells (Figures SID and S1E). After 
several days of tissue culture, we measured and compared 
the expression profiles of these engineered tissues. A com- 
parison of the expression profiles for anisotropic neonate 
and mESC tissues (Figure IDi) revealed a number of differ- 
ences associated with Ca 2+ channel subunits, such as the 
L-type Ca 2+ channel subunit Cacnald (37.5-fold, p < 
0.0001), as well as the T-type subunits Cacnalg (20.2-fold, 
p < 0.05) and Cacnalh (23.8-fold, p < 0.05). Additionally, 
the mESC anisotropic tissues showed significantly lower 
expression of the ventricular marker Irx4 (—7.7-fold, 
p < 0.05) and significantly higher expression of the atrial 
markers Myl4 (254.8-fold, p < 0.01) and Myl7 (104.0-fold, 
p < 0.01) versus the neonate tissues. 

In contrast, the miPSC anisotropic tissues exhibited sig- 
nificant differences from the neonate tissues (Figure IDii) 
for the Ca 2+ channel subunits Cacnald (36.9-fold, p < 
0.05) and Cacnalg (6.6-fold, p < 0.05), as well as the atrial 
myosin light chain kinase gene Myl4 (105.5-fold, p < 
0.01). Hierarchical clustering of neonate, mESC, and 
miPSC gene-expression measurements revealed a distinct 
separation of the expression profiles for isotropic and 
anisotropic tissues, regardless of myocyte type (Figure IE). 
Moreover, the expression profiles for mESC and miPSC 
myocytes in both the isotropic and anisotropic cellular 
configurations clustered closer to each other than to the 
neonate tissues. This suggests that the mESC and miPSC 
myocytes exhibited global transcriptional profiles that 
were distinct from the neonate expression pattern, despite 
differences in the relative expression profiles between the 
mESC and miPSC tissues. 

Characterization of Myofibril Architecture and Global 
Sarcomere Alignment 

One of the defining features of the ventricular myocardium 
is the laminar arrangement of cardiac myocytes, which 
serves to organize and orient the contractile sarcomeres 
to facilitate efficient pump function (McCain and Parker, 
2011). We evaluated the ability of the mESC and miPSC 
engineered tissues to self-assemble myofibrils with align- 
ment comparable to that of neonate ventricular myocytes 
by using custom image analysis software developed in 
our lab. Immunofluorescence micrographs of sarcomeric 
a-actinin allowed us to visualize the orientations of the z 
lines outlining the lateral edges of sarcomeres, and to quan- 
titatively assess sarcomere organization in our engineered 
tissues (Figure S2Bi-S2Bii). Visualization of global z-line 
registration in isotropic monolayers of mESC (Figure 2Ai), 
miPSC (Figure 2Aii), and neonate (Figure 2Aiii) myocytes 
revealed random orientation patterns. In contrast, the 
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Figure 1. Comparison of mESC, miPSC, and Neonate Gene-Expression Profiles on Isotropic and Anisotropic ECM Substrates 

(A) CuLturing (i) mESC, (ii) miPSC, and (iii) neonate myocytes on substrates with a uniform coating of FN resulted in isotropic cellular 
arrangement. 

(B) Volcano plots showing the negative log of p values (two-tailed t test, n = 3 for all conditions) versus log fold-change values for 
comparison of qPCR measurements of cardiac genes (i) between mESC and neonate isotropic monolayers, and (ii) between miPSC and 
neonate isotropic monolayers reveal significant differences for a number of genes (points on the plot colored green or red represent genes 
with p < 0.05). 

(C) Culturing (i) mESC, (ii) miPSC, and (iii) neonate myocytes on substrates with microcontact-printed lines of FN that were 20 urn wide 
and spaced 4 urn apart resulted in anisotropic cellular arrangement in all three cell types. 

(D) Volcano plots showing the negative log of p values (two-tailed t test, n = 3 for aLl conditions) versus log fold-change values for 
comparison of qPCR measurements of cardiac genes (i) between mESC and neonate anisotropic monolayers, and (ii) between miPSC and 
neonate anisotropic monolayers reveal slightly fewer genes demonstrating significant differences than in the isotropic cultures (points on 
the plot colored green or red represent genes with p < 0.05). 

(legend continued on next page) 
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anisotropic mESC (Figure 2Bi), miPSC (Figure 2Bii), and 
neonate (Figure 2Biii) tissues demonstrated a greater degree 
of uniaxial z-line registration. To quantify the differences in 
global sarcomere organization between the mESC and 
miPSC tissues versus the neonate tissues (Figure 2C), we 
utilized a metric known as the orientational order param- 
eter (OOP), which is commonly used to characterize the 
alignment of liquid crystals (Kuczynski et al., 2002) and 
ranges from zero (random organization) to one (perfect 
alignment). It was previously shown that this metric can 
be successfully adapted to measure and compare z-line 
registration in engineered cardiac tissues, and that it pro- 
vides insight into the contractile strength of the tissues 
(Feinberg et al., 2012). The anisotropic neonate tissues ex- 
hibited a significantly higher OOP value than both the 
mESC and miPSC tissues, suggesting that both types of 
stem cell-derived cardiac myocytes were unable to generate 
myofibrils with the same degree of global sarcomere align- 
ment as the neonate myocytes. Isotropic tissues had low 
OOP values due to the random organization of the cardiac 
myocytes. Measurement of registered z-line spacing also re- 
vealed that the anisotropic mESC and miPSC tissues dis- 
played significantly shorter sarcomere lengths than the 
neonate tissues (Figure 2D). Moreover, quantification of 
"sarcomere packing density" (i.e., the proportion of a-acti- 
nin localized to z lines indicative of the presence of fully 
formed sarcomeres) showed that the anisotropic neonate 
tissues exhibited significantly higher sarcomere packing 
density than the mESC and miPSC tissues. Taken together, 
these analyses revealed that the mESC- and miPSC-derived 
myocytes responded to ECM cues in a manner similar to 
that observed for the neonate myocytes, but exhibited 
sarcomere organization reminiscent of immature premyo- 
fibrils observed in embryonic cardiac myocytes (Dabiri 
et al, 1997; Grosberg et al., 2011b; LoRusso et al., 1997). 

Measurement of mESC, miPSC, and Neonate 
Electrophysiological Performance 

The electrical activity of cardiac myocytes regulates the 
initiation of myofibril contraction and is commonly 
measured as an indicator of myocyte identity and function- 
ality (Kleber and Rudy, 2004; Maltsev et al., 1994; Weinberg 
et al., 2010). We used planar patch-clamp recordings to 
compare and contrast the action potential (AP) characteris- 
tics of isolated mESC, miPSC, and neonate myocytes. We 
identified two different demographics of cell types as 
demonstrated by AP morphology (Maltsev et al., 1994). 



Neonate myocytes primarily demonstrated ventricular- 
like APs (Figure 3Ai), whereas mESC- and miPSC-derived 
myocytes exhibited APs that were evenly distributed be- 
tween ventricular-like and atrial-like morphologies (Fig- 
ure 3Aii). Both the mESC- and miPSC-derived myocytes 
primarily exhibited APs, as shown in Figure 3Aii, whereas 
the neonate ventricular myocytes demonstrated APs, as 
illustrated in Figure 3Ai. Analysis of AP characteristics, 
such as maximum voltage (V max ), AP duration at 50% repo- 
larization (APD50), and AP duration at 90% repolarization 
(APD90), revealed that the mESC and miPSC myocytes 
exhibited roughly equal incidences of atrial-like and ven- 
tricular-like APs, whereas the neonate cardiac myocytes dis- 
played ventricular-like AP characteristics (Figure 3B). In 
addition to AP characterization, we also measured the elec- 
trical conduction properties of the anisotropic mESC (Fig- 
ure 3Ci), miPSC (Figure 3Cii), and neonate (Figure 3Ciii) 
tissues using optical mapping and the voltage-sensitive 
fluorescent dye RH-237 (Bursac et al., 2002; Thomas 
et al., 2000; Weinberg et al., 2010) to evaluate the ability 
of the stem cell-derived myocytes to form the electrome- 
chanical syncytium that typifies the myocardium (Kleber 
and Rudy, 2004). We did not observe any significant differ- 
ences in the longitudinal (LCV) or transverse (TCV) 
conduction velocities between the mESC, miPSC, and 
neonate tissues (Figure 3D). However, we did observe sub- 
stantial differences in the cellular dimensions (/ x w) of 
mESC (72.04 um x 12.07 um) and miPSC (82.30 urn x 
11.25 urn) versus neonate (44.93 um x 11.35 um) cardiac 
myocytes that may influence the magnitude of the LCV 
and TCV in engineered tissues comprised of these 
cells, making them appear faster than they actually are. 
AP duration measurements revealed no significant differ- 
ences at 50% repolarization (APD50), but a significant 
(p < 0.05) difference was observed at 90% repolarization 
(APD90) between the neonate and mESC anisotropic tis- 
sues (Figure 3E). 

Ca 2+ plays a crucial role in coupling myocyte excitation 
and contractile activity (Bers, 2002). We measured Ca 2+ 
transient activity in engineered anisotropic tissues, as 
well as the Ca 2+ current profiles of isolated mESC, miPSC, 
and neonate myocytes. Ca 2+ transients measured in aniso- 
tropic tissues revealed a significantly (p < 0.05) shorter 50% 
decay time in the miPSC, but not the mESC, tissues as 
compared with the neonate tissues, and significantly 
(p < 0.05) shorter 90% decay time in both the mESC and 
miPSC tissues versus the neonate tissues (Figure 3F). Planar 



(E) Hierarchical clustering of mean 2 values for a select panel of cardiac genes reveals that the isotropic and anisotropic neonate tissue 
expression profiLes duster together in the center columns of the heatmap, whereas the anisotropic mESC and miPSC expression profiles 
form a separate cluster on the right sides of the heatmap, and the isotropic mESC and miPSC profiles cluster together on the left side of the 
heatmap. Scale bars, 100 um. 
See also Figure SI and Table SI. 
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Figure 2. Comparison of Myofibril Architecture in mESC, miPSC, and Neonate Engineered Tissues 

(A and B) Immunofluorescence visualization of sarcomeric a-actinin in (A) isotropic monolayers of (i) mESC, (ii) miPSC, and (iii) neonate 
myocytes, and (B) anisotropic monolayers of (i) mESC, (ii) miPSC, and (iii) neonate myocytes reveals the pattern of sarcomere organization 
adopted by each cell type in response to geometric cues encoded in the ECM. Immature premyofibrils (red arrows) were observed 
exclusively in mESC and miPSC engineered tissues. Quantitative evaluation of sarcomeric a-actinin immunofluorescence micrographs 
allowed statistical comparison of sarcomere organization and architecture. 

(C) The OOP was used as a metric of global sarcomere alignment within the engineered tissues and showed that anisotropic neonate tissues 
exhibited significantLy greater overall sarcomere alignment than the mESC and miPSC anisotropic tissues. No significant differences in 
global sarcomere alignment were observed among the isotropic mESC, miPSC, and neonate tissues. 

(D) Comparison of z-line spacing revealed that the neonate anisotropic tissues exhibited significantly greater sarcomere length than both 
the mESC and miPSC anisotropic tissues. 

(E) From the measurements of sarcomere length, the sarcomere packing density was calculated for anisotropic tissues of each cell type. All 
three cell types exhibited significantly different sarcomere packing densities. The statistical tests used were ANOVA (*p <0.05) and ANOVA 
on ranks (+p < 0.05). Data are presented as mean + SEM. ScaLe bars, 10 u.m. 

See also Figure S2. 

patch-clamp recordings of L-type (LCC; Figure 3Gi) and versus the mESC-derived, but not the miPSC-derived, myo- 
T-type (TCC; Figure 3Gii) Ca 2+ current profiles revealed cytes (Figure 3H). Taken together, these data suggest that 
significantly (p < 0.05) higher total (TOT) and TCC the mESC and miPSC myocytes possessed electrophysio- 
maximum Ca 2+ current densities in the neonate myocytes logical properties similar to those of neonate cardiac 
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myocytes, aside from differences in funny current and 
voltage-gated Ca 2+ channel subunit expression (illustrated 
in Figure 1). 

Measurement of mESC, miPSC, and Neonate 
Engineered Tissue Contractile Performance 
With the muscular thin film (MTF) contractility assay it is 
now possible to assess the diastolic (Figure 4Ai) and systolic 
(Figure 4Aii) functions of engineered myocardium directly 
(Alford et al., 2010; Feinberg et al. ( 2007; Grosberg et al., 
2011a). Using the "heart-on-a-chip" MTF assay (Grosberg 
et al., 2011a), we measured the stress generation profiles 
of anisotropic mESC, miPSC, and neonate tissues (Fig- 
ure 4B), and compared their contractile performance. The 
anisotropic neonate tissues generated significantly (p < 
0.05) higher diastolic, peak systolic, and twitch stress 
than both the mESC and miPSC tissues (Figure 4C), with 
the observed values for the neonate tissues being within 
the range measured for isolated murine papillary muscle 
strips (Gao et al., 1998; Stuyvers et al., 2002). The results 
of the contractility measurements clearly show a functional 
deficit in the mESC- and miPSC-derived myocytes that was 
not apparent in the electrophysiological measurements. 
We used the combined output of our electrophysiological, 
calcium transient, and contractile force experimental mea- 
surements to create graphical representations of the excita- 
tion-contraction coupling profiles of the mESC (Figure 4Di), 
miPSC (Figure 4Dii), and neonate (Figure 4Diii) engineered 
tissues that clearly illustrate the similarities and differences 
in the excitation-contraction coupling among the cell 
types. These data show that the miPSC-derived myocytes 
are qualitatively more similar to the neonate myocytes 
than are the mESC-derived myocytes. 

Integration of Experimental Measurements to 
Evaluate the Stem Cell-Derived Myocyte Cardiac 
Phenotype 

To determine how closely the mESC- and miPSC-derived 
myocytes matched the phenotype of the neonate ventricu- 
lar myocytes, we computationally integrated the set of 
gene-expression, morphology, electrophysiology, and con- 
tractility experimental measurements collected from each 
cell population, and calculated the difference between 
the unknown and target cell populations. For this purpose, 
we evaluated measures of effect size, such as the z factor 
and strictly standardized mean difference (SSMD), that 
have previously been used to quantify biological popula- 
tion differences in high-throughput screening applications 
(Birmingham et al., 2009). We chose to use the SSMD 
instead of the z factor in this case study because the 
SSMD is more robust to outliers, is not dependent on sam- 
ple size, and can be used to evaluate non-normal data 
(Zhang, 2007). For each type of experimental measure- 



ment, we normalized the mean (iinorm) values to the inter- 
val [0, 1] by identifying the maximum (u. max ) and minimum 
(Umin) mean values observed for that measurement from all 
three cardiac myocyte sources, and calculating as follows: 

M„ am ,= (Equation 1) 



Using these normalized values, we calculated the SSMD 
(P) and quantified the differences between each unknown 
population (i.e., mESC and miPSC) and the neonate target 
population as follows: 



^target ^unknown 



ij2 I _-2 

target ^ unknown 



(Equation 2) 



where u. represents the mean and a represents the SD. We 
used the results to evaluate the magnitude of difference, 
taking into account the variance in the measurements, 
between the stem cell-derived myocytes and the neonate 
cardiac myocytes (Figure 5). This allowed us to identify 
the parameters that showed the greatest degree of similarity 
to and difference from the target neonate ventricular myo- 
cyte tissues. We then used the p values from each experi- 
mental measurement for the mESC and miPSC tissues, 
and calculated the mean squared error (MSE) versus the 
neonate tissues as follows: 



MSE-- 



1 



n 4—? 



(Equation 3) 



where n is the total number of experimental measurement 
P values included in the calculation. We then used these 
results to evaluate the differences observed for each mea- 
surement category (i.e., the p values for gene expression, 
morphology, electrical activity, contractility used to calcu- 
late category-specific MSE values), as well as to define a 
single MSE value calculated from all of the experimental 
measurements from all categories combined that repre- 
sents the total difference between the stem cell-derived 
and neonate cardiac myocytes based on the measurements 
performed (Table 1). A lower MSE value indicates a better 
match to the neonate target phenotype, with an MSE value 
of zero indicating a perfect match. 

We found that the miPSC tissues exhibited lower MSE 
values than the mESC tissues in every measurement cate- 
gory except morphology. In addition, the overall MSE 
values calculated from all of the experimental measure- 
ments combined revealed a lower MSE for the miPSC engi- 
neered tissues than for those comprised of mESC-derived 
myocytes. This suggests that the miPSC-derived myocytes 
exhibited a global phenotype that was slightly closer to 
the neonate cardiac myocytes than the mESC-derived 
myocytes, although both the mESC- and miPSC-derived 
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Figure 3. Comparison of Electrical Activity in mESC, miPSC, and Neonate Engineered Tissues 

(A) Patch-clamp recordings from isolated mESC, miPSC, and neonate myocytes exhibited APs with both (i) ventricular-like and (ii) atrial- 
like profiles. 

(B) Characterization of the AP traces revealed no significant differences between the three cell types, but the mESC and miPSC myocytes 
exhibited an equal proportion of ventricuLar-like (mESC-v, miPSC-v) and atrial-like (mESC-a, miPSC-a) AP traces, whereas the neonates 
exhibited primarily ventricular-like (neonate-v) AP profiles. 

(C) The electrophysiological characteristics of anisotropic (i) mESC, (ii) miPSC, and (iii) neonate tissues were assessed using optical 
mapping and the photovoltaic dye RH237. 

(D) Comparison of conduction properties between the mESC, miPSC, and neonate tissues revealed no significant differences in either 
LCV or TCV. 

(E) Evaluation of optical AP duration in anisotropic tissues revealed no significant differences in APD50, but a significant difference in 
APD90 between the mESC and neonate tissues was observed. 

(legend continued on next page) 
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Figure 4. Comparison of Contractile Per- 
formance in mESC, miPSC, and Neonate 
Engineered Tissues 

(A) The contractile performance of aniso- 
tropic mESC, miPSC, and neonate tissues 
was assessed using the MTF assay, and the 
radius of curvature of the MTFs at (i) dias- 
tole and (ii) peak systole were used to 
calculate contractile stress. 

(B) The radius of curvature of the MTFs was 
used to calculate and compare the temporal 
contractile strength profiles of anisotropic 
mESC (green), miPSC (red), and neonate 
(blue) tissues. 

(C) Comparison of MTF contractile output 
revealed that neonate anisotropic tissues 
generated significantly greater diastolic, 
peak systolic, and twitch stress than both 
the mESC and miPSC tissues. 

(D) Graphical representation of AP mor- 
phology (black solid line), Ca 2+ transient 
morphology (blue dotted line), and con- 
tractility profile (red dotted line) during a 
typical excitation-contraction cycle of the 
mESC, miPSC, and neonate engineered 
anisotropic tissues. The statistical test used 
was AN0VA (*p < 0.05). Data are presented 
as mean ± SEM. 

See also Movies SI, S2, and S3. 
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myocytes demonstrated substantial differences from the 
neonate cardiac myocytes for a number of characteristics. 

DISCUSSION 

The goal of this study was to define a quality-control stan- 
dard rubric for assessing stem cell-derived cardiac myo- 
cytes. We chose a set of experimental measurements that 
provide insight into not only the expression profile of the 



cells but also morphological and functional characteristics 
that are intimately tied to the contractile function of car- 
diac tissues (Bursac et al., 2002; Feinberg et al. ( 2012; Kleber 
and Rudy 2004). We utilized ventricular myocytes isolated 
from postnatal mouse hearts to serve as our reference stan- 
dard for denning the target phenotype. However, an 
inherent limitation to using these cells is the presence of 
noncardiomyocytes, such as fibroblasts, endothelial cells, 
and smooth muscle cells, that may confound the inter- 
pretation of some experimental measurements, such as 



(F) Comparison of Ca 2+ transients measured in anisotropic tissues revealed that the 50% decay time of the miPSC tissues was significantly 
lower than that of both the mESC and neonate tissues, but the 90% decay time of both the mESC and miPSC tissues was significantly lower 
than that of the neonate tissues. 

(G) Patch-clamp recordings were collected on isolated mESC, miPSC, and neonate myocytes to measure and compare (i) LCC and (ii) TCC 
elicited at various holding potentials. 

(H) Patch-clamp recordings of maximum Ca 2+ current density in isolated mESC, miPSC, and neonate myocytes revealed a significant 
difference in TOT between the neonate and mESC myocytes. No significant differences in LCC were observed, but a significant difference 
in TCC was observed between the neonate and mESC myocytes. The statistical test used was AN0VA (*p < 0.05). Data are presented as 
mean ± SEM. Scale bars, 20 urn. 

See also Figure S3. 
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Figure 5. Integrated Visual Comparison of mESC, miPSC, and 
Neonate Experimental Measurements 

SSMD (P) vaLues were computed for mESC- and miPSC-derived 
myocytes relative to the neonate cardiac myocytes from the mean 
and sample SDs collected for each experimental measurement. 
These (3 values were organized by measurement type (i.e., gene 
expression, myocyte architecture, electrophysiology, and con- 
tractility) and plotted to allow comparison. Negative fi values 
indicate measurements with higher relative magnitude in the 
neonate cardiac myocytes, whereas positive fi values indicate 
measurements that were higher in the mES/miPSC myocytes rela- 
tive to the neonate cardiac myocytes. 
See also Table S2. 



gene-expression profiling. It should also be noted that the 
mESC- and miPSC-derived cardiac myocytes used in this 
study were produced using a differentiation protocol that 
gives rise to a heterogeneous population of atrial, ventricu- 
lar, and pacemaker-like cells, which may also make inter- 
pretation of certain measurements challenging. However, 
the proposed quality-assessment strategy is not dependent 
on any particular set of measurements and allows re- 
searchers the flexibility to choose the set of experimental 
measurements that best suits their needs. 

Using the experimental measurements described above 
and isolated neonatal ventricular myocytes as our reference 
phenotype, we developed a "quality index" that utilizes the 



Table 1. MSE Values Calculated for Each Group of Measurements 
in the Comparison of mESC- and miPSC-Derived Myocytes with 
Neonate Ventricular Myocytes 


Measurement Category 


MSE mES 


MSE miPS 


Gene expression 


5.69 


4.25 


Morphology 


1.30 


1.48 


Electrophysiology 


1.16 


0.57 


Contractility 


6.32 


2.95 


All measurements 


4.95 


3.60 



The SSMD (P) values computed for each experimental measurement were 
used to calculate MSE values for each major measurement category (see 
Table S2), as well as all of the measurements combined, in comparisons of 
the mESC (MSE mES ), and miPSC (MSE m j PS ) engineered tissues with the 
neonate engineered tissues. 



magnitude and variance of these measurements to provide 
a numeric score indicating how closely the stem cell- 
derived myocytes match the characteristics of the neonatal 
cardiac myocytes. The combination of gene-expression, 
morphological, electrophysiological, and contractility 
measurements employed allowed us to pinpoint specific 
differences in the structural and functional properties of 
the mESC and miPSC engineered tissues versus the neonate 
tissues that have important implications for their utility in 
in vitro assays. Additional studies of the relationships be- 
tween these measurements and the response of engineered 
cardiac tissues to compounds that have known effects on 
heart function may provide valuable insight into the com- 
bination of measurements that can most reliably determine 
the ability of stem cell-derived cardiac myocytes to adopt 
the desired phenotype. With a carefully chosen set of exper- 
imental parameters, this quality-assessment rubric may 
provide a reliable means to evaluate strategies for improving 
the differentiation of cardiac myocytes from stem cells and 
drive them toward a more mature phenotype (Cassino et al., 
2012). Further, this quality index will not only allow re- 
searchers to identify the commercial stem cell-derived 
myocyte product lines that are most suitable for their needs, 
it may also serve the stem cell industry as a quality-assur- 
ance system for ensuring that batches released to customers 
faithfully recapitulate the desired phenotype. 

EXPERIMENTAL PROCEDURES 

Stem Cell-Derived Myocyte Culture 

Cor.At mESC-derived myocytes (lot* CS25CL_V_SN_1M, produc- 
tion date: 2010-02-16) and iPSC-derived myocytes (lot* CS02CL-i, 
production date: 2010-11-15; CS07CL4, production date: 2011-05- 
09) were cultured according to instructions and with culture 
reagents supplied by the manufacturer (Axiogenesis). Briefly, 
cells were cultured in T25 flasks precoated with 0.01 ng/ml FN 
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(BD Biosciences) in puromycin-containing culture media at 37°C 
and 5% C0 2 for 72 hr to eliminate undifferentiated stem cells 
from the culture. After 72 hr, cells were dissociated with 0.25% 
trypsin and seeded onto microcontact-printed substrates at 
densities of 100,000/cm 2 . Cells were cultured for 2 days on micro- 
contact-printed substrates to allow formation of a functional 
syncytium while preserving the proportion of myocytes that 
exhibited ventricular characteristics prior to experimentation. 

Neonatal Mouse Ventricular Myocyte Culture 

Neonatal mouse ventricular myocytes were isolated from 2-day-old 
neonatal Balb/c mice using procedures approved by the Harvard 
University Animal Care and Use Committee. Briefly, excised ven- 
tricular tissue was incubated in a 0.1% (w/v) trypsin (USB) solution 
cooled to 4°C for approximately 12 hr with agitation. Trypsinized 
ventricular tissue was dissociated into cellular constituents via 
serial exposure to a 0. 1% (w/v) solution of collagenase type II (Wor- 
thington Biochemical) at 37°C for 2 min. Isolated myocytes were 
maintained in a culture medium consisting of Medium 199 (Invi- 
trogen) supplemented with 10% (v/v) heat-inactivated fetal bovine 
serum (FBS), 10 mM HEPES, 20 mM glucose, 2 mM L-glutamine, 
1.5 uM vitamin B-12, and 50 U/ml penicillin, and seeded at a den- 
sity of 200,000 cells/cm 2 . From the second day of culture onward, 
the FBS concentration was reduced to 2% (v/v) and the medium 
was exchanged every 48 hr. Myocytes were cultured for 4 days on 
microcontact-printed substrates prior to experimentation. 

Fabrication of Microcontact-Printed Substrates 

Silicone stamps designed for microcontact printing were prepared 
as previously described. Photolithographic masks were designed in 
AutoCAD (Autodesk), and consisted of 20 urn wide lines separated 
by 4 (im gaps to impose a laminar organization on the myocytes. 
Polydimethylsiloxane (PDMS, Sylgard 184; Dow Corning) was 
used to fabricate stamps with the specified pattern. Stamps were 
incubated with 50 ng/ml FN (BD Biosciences) for 1 hr. Glass cover- 
slips were spin-coated with PDMS and treated in a UV-ozone 
cleaner (Jelight) immediately prior to stamping with FN. After 
the FN pattern was transferred to the surface of the PDMS-coated 
coverslips, they were incubated in 1% (w/v) Pluronic F-127 
(BASF) to block cell adhesion to unstamped regions. 

Heart-on-a-Chip Substrate Fabrication 

Engineered cardiac tissue contractile performance was measured 
using a custom MTF based platform (Alford et al., 2010; Feinberg 
et al., 2007; Grosberg et al., 2011a). Briefly, the heart-on-a-chip 
substrates consisted of glass coverslips selectively coated with a 
thermosensitive sacrificial polymer, poly(N-isopropylacrylamide) 
(PiPAAm; Polysciences), and with a second layer of PDMS. The 
thickness of the PDMS layer was found to be in the range of 10- 
18 urn for all "heart chips" used in this study (Dektak 6M; Veeco 
Instruments). 

Heart-on-a-Chip Contractility Experiments 

During contractility experiments, samples were submerged in 
Tyrode's solution (mM, 5.0 HEPES, 5.0 glucose, 1.8 CaCl 2 , 1.0 
MgCl 2 , 5.4 KC1, 135.0 NaCl, and 0.33 NaH 2 P0 4 , pH 7.4). All re- 
agents were purchased from Sigma-Aldrich. Rectangular films 



were cut with a razor blade, and the bath temperature was 
decreased below the PiPAAm transition temperature, making 
possible for the MTF to bend away from the glass. Video recordings 
of the deformation of each film were processed to obtain the time 
course (Alford et al., 2010) of the tissue-generated stresses. The 
peak systolic and diastolic stresses were calculated as the average 
of the maxima and minima of the stress profile during 10 cycles 
at a pacing of 3 Hz, and twitch stress was defined as the difference 
between peak systolic and diastolic stresses. 

Immunohistochemical Labeling 

Samples were fixed in 4% (v/v) paraformaldehyde with 0.05% (v/v) 
Triton X-100 in PBS at room temperature for 10 min. Cells were 
incubated in a solution containing 1:200 dilutions of monoclonal 
anti-sarcomeric a-actinin antibody (A7811, clone EA-53; Sigma- 
Aldrich), polyclonal anti-FN antibody (F3648; Sigma-Aldrich), 
DAPI (Invitrogen), and Alexa Fluor 633-conjugated phalloidin 
(Invitrogen) for 1 hr at room temperature. Samples were then incu- 
bated in 1:200 dilutions of Alexa Fluor 488-conjugated goat anti- 
mouse IgG and Alexa Fluor 546-conjugated goat anti-rabbit IgG 
secondary antibodies (Invitrogen) for 1 hr at room temperature. 
Labeled samples were imaged with a Zeiss LSM confocal micro- 
scope (Carl Zeiss Microscopy). 

Quantitative Evaluation of Sarcomere Structure 

Confocal Z stacks of sarcomeric a-actinin fluorescence micrographs 
were deconvolved with Mediacy Autoquant (MediaCybernetics). 
Analysis of sarcomeric structural characteristics was conducted 
using custom-designed ImageJ (Abramoff et al., 2004) and MATLAB 
(The Math Works) software (Figure S2A). The orientations of sarco- 
meric a-actinin positive pixels were determined using an adapted 
structure-tensor method (Rezakhaniha et al., 2012) and the orienta- 
tional order parameter (OOP), a measure of the global alignment of 
the sarcomeres, was calculated from the observed orientation 
values. The orientations observed in the micrographs were color- 
coded using the HSV digital image representation (Figure S2Bi) 
where the Hue channel was used for orientation, the Saturation 
channel for pixel coherency (i.e., a measure of local contrast), and 
the Value channel for the preprocessed image. The normalized 
occurrence of orientations that demonstrated a coherency higher 
than a given threshold could then be displayed in a histogram (Fig- 
ure S2Bii). The sarcomere length and the overall regularity of the 
z lines were determined by processing the images with a 2D fast 
Fourier transform algorithm. The power spectrum (Figure S2Ci) 
was then radially integrated and normalized by the total area under 
the ID curve. The sarcomere packing density was defined as the 
area under the periodic component (shaded red in Figure S2Cii). 

Planar Patch-Clamp Electrophysiological Recordings 

Planar patch-clamp experiments were conducted as previously 
described (Weinberg et al., 2010). Briefly, cells were cultured on 
FN (BD Biosciences) coated T25 flasks for 5 days, isolated using 
0.25% trypsin (Invitrogen), resuspended in extracellular buffer so- 
lution (EBS [mM]: 140 NaCl, 4 KC1, 1 MgCl 2 , 2 CaCl 2 , 5 D-glucose 
monohydrate, 10 HEPES, pH 7.4) to a final concentration of 1,000 
cells/ ul, and allowed to equilibrate for 5 min in EBS. The electronics 
were calibrated in the presence of EBS and intracellular buffer 
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solution (IBS [mM] : 50 KC1, 10 NaCl, 60 KF, 20 EGTA, 10 HEPES, pH 
7.2) prior to flowing cells into the chamber. Then 5 ul of the cell sus- 
pension was introduced into the chip and the negative pressure 
was automatically adjusted to produce a final seal resistance greater 
than 1 GOhm. During current-clamp experiments, cells were sub- 
jected to ten trains of ten current pulses at 3 Hz; the current ampli- 
tude was set to 1 . 5 times the threshold for AP generation. When the 
signal reached steady state, ten APs were averaged, yielding a repre- 
sentative trace for the calculation of AP duration indicators. During 
voltage-clamp experiments, cells were kept in buffers containing 
tetrodotoxin (10 uM), nifedipine (10 jiM), 4-AP (1 mM), and tetrae- 
thylammonium (20 mM) purchased from Sigma- Aldrich. The 
membrane potential was subjected to two voltage-clamp protocols: 
First, the membrane potential was held to a value of -90 V for 
250 ms and then stepped from -70 to +40 mV in 10 mV steps for 
250 ms, thus eliciting the total Ca 2+ current (TOT). Second, from 
the same holding potential, cells were stepped from -40 
to +40 mV, a range in which mostly the LCC is active. The TCC 
was then calculated as the difference between the TOT and LCC. 

Optical Mapping of Electrophysiological Properties 

Samples were incubated in 4 uM RH237 (Invitrogen) for 5 min and 
washed three times with Tyrode's solution prior to recording. The 
temperature of the bath solution was maintained at approximately 
35°C using a digital temperature controller (TC-344B; Warner In- 
struments) for the duration of the experiment. Then 10 uM Bleb- 
bistatin (EMD Millipore) was added to minimize motion artifacts 
during recording of electrical activity (Fedorov et al., 2007; Herron 
et al., 2012). Samples were paced at 3 Hz with a 10 ms biphasic 
pulse at 10-15 V delivered using an SD-9 stimulator (Grass Tech- 
nologies) and a bipolar, platinum point electrode placed approxi- 
mately 300-500 jim above the sample and 1-2 mm from the 
top-right corner of the field of view (FOV). Imaging was performed 
using a Zeiss Axiovert 200 epifluorescence microscope (Carl Zeiss 
Microscopy) equipped with an X-cite Exacte mercury arc lamp 
(Lumen Dynamics). Illumination light was passed through a 
40X/1.3 NA objective (EC Plan-NEOFLUAR; Zeiss) and a band- 
pass excitation filter (530-585 nm). Emission light was filtered at 
615 nm with a long-pass filter and focused onto the 100 x 100 
pixel chip of a high-speed MiCAM Ultima CMOS camera (Scime- 
dia). Images were acquired at 1,000 fps from a 250 x 250 nm 
FOV. Postprocessing of the raw data included reduction of drift 
induced by photobleaching by subtracting a linear fit of the base- 
line, applying a 3 x 3 pixel spatial filter to improve the signal/noise 
ratio, and excluding saturated pixels. Activation time was calcu- 
lated as the average maximum upstroke slope of multiple pulses 
over a 2-4 s recording window. LCVs and TCVs were calculated 
through a linear fit of the activation times along the horizontal 
and vertical axes of each FOV, respectively. Optical AP traces 
were calculated as the average of multiple pulses while adjusting 
the offset of each pixel caused by different activation times. 

Ratiometric Measurement of Cardiac Myocyte Ca 2+ 
Transients 

Cardiac tissues were incubated in a 5 uM solution of acetoxy- 
methyl (AM) Fura Red (F-3021; Invitrogen) reconstituted in 
Pluronic F-127 (P-3000MP; Invitrogen) for 20 min. After dye 



loading, cells were incubated in Tyrode's solution for 5 min, rinsed 
three times, and imaged on a Zeiss LSM LIVE (Carl Zeiss Micro- 
scopy) confocal microscope at 40 x magnification with an environ- 
mental chamber to ensure a constant physiological temperature of 
37°C. Tissues were field stimulated at 3 Hz during recordings. Dual- 
excitation ratiometric recordings were performed by rapidly 
switching (through an acousto-optical tunable filter) excitation 
laser lights at 405 nm and 488 nm, and collecting the correspond- 
ing emissions through a high-pass filter with cutoff at 546 nm. 
Recordings were constrained to 20 lines oriented perpendicular 
to the main axis of the cells and ensuring minimal intersection 
with nuclei to maintain an acquisition speed of 250 fps (Fig- 
ure S3 A). After background subtraction, two signals were obtained 
(Figure S3B): one (blue line) that increased with the Ca 2+ elevation 
corresponding to excitation at 405 nm, and one (green line) that 
showed an opposite trend and corresponded to the 488 nm excita- 
tion wavelength. The ratiometric representation of the Ca 2+ tran- 
sient was taken as the ratio of the 405 nm and 488 nm signals 
(black trace in Figure S3C). Four consecutive transients at steady 
state were further averaged to create a representative single tran- 
sient (Figure S3D). 

Quantitative RT-PCR Gene-Expression Measurements 

Total RNA was collected in triplicate from both isotropic and 
micropatterned anisotropic samples using the Strategene Abso- 
lutely RNA Miniprep kit (Agilent Technologies) according to the 
manufacturer's instructions. Genomic DNA contamination was 
eliminated by incubating the RNA lysates in DNase I digestion 
buffer at 37°C for 15 min during the RNA purification procedure. 
The quantity and purity of the RNA lysates were assessed using a 
Nanodrop spectrophotometer (Thermo Scientific). Purified total 
RNA lysates with OD 260/280 ratios greater than 1.7 were used 
for quantitative RT-PCR (qRT-PCR) measurements. cDNA strands 
were synthesized for genes of interest using an RT2 first-strand syn- 
thesis kit (QIAGEN) and custom preamplification primer sets 
(QIAGEN). For each first-strand synthesis reaction, 500 ng of total 
RNA from each lysate was used. Expression levels for specific genes 
of interest (Table SI) were measured using custom RT2 Profiler RT- 
PCR arrays (QIAGEN) and the Bio-Rad CFX96 RT-PCR detection 
system. Statistical analysis of qRT-PCR threshold cycle data was 
carried out with the web-based RT2 Profiler PCR Array Data anal- 
ysis Suite version 3.5 (QIAGEN) according to published guidelines 
(Schmittgen and Livak, 2008). 

SUPPLEMENTAL INFORMATION 

Supplemental Information includes Supplemental Experimental 
Procedures, three figures, two tables, and three movies and can 
be found with this article online at http://dx.doi.Org/10.1016/j. 
stemcr.2014.01.015. 
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